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f**^ . We apply the effective potential analytic continuation (EPAC) method to one-dimensional asym- 

metric potential systems to obtain the real time quantum correlation functions at various tem- 
peratures. Comparing the EPAC results with the exact results, we find that for an asymmetric 



anharmonic oscillator the EPAC results are in very good agreement with the exact ones at low 
temperature, while this agreement becomes worse as the temperature increases. We also show that 
the EPAC calculation for a certain type of asymmetric potentials can be reduced to that for the 
corresponding symmetric potentials. 
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I. INTRODUCTION 
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The imaginary time path integral formalism [1] has been an important tool for studying static or dynamical 
properties of quantum statistical mechanical systems. The static properties have successfully been calculated by means 
of the path integral Monte Carlo (PIMC) or path integral molecular dynamics (PIMD) technique HQ- However, these 



numerical methods have difficulties in computing the dynamical properties such as the real time quantum correlation 
functions, because the analytic continuation from the imaginary time to the real time is a non-trivial procedure in 



the numerical calculations 0, 

As new quantum dynamics methods to calculate the real time quantum correlation functions at finite temperature, 
the centroid molecular dynamics (CMD) method Q and the effective potential analytic continuation (EPAC) method 
have recently been proposed. Both the methods use effective potentials, which can be numerically calculated by 
means of the PIMC or PIMD technique. In computing the real time quantum correlation functions, the CMD method 

lfj, while the EPAC method uses the standard effective potential which 

tials have been extensively used as powerful 



uses the effective classical potential 



111 ll2jL The effective pot entia 
!Tties|a E3 EJ El El Q , so 



appears in the effective action formalism 

tools in the calculations of the static properties [9|, llCj, |llj, 1 121 1 1.1 Il4| , so that they are expected to be effective also 
for the calculations of the dynamical properties. In practice, the CMD method has been tested in low-dimensional 
systems H 0, 0, Q and applied to various many-body molecular systems to yield many promising results 
|l9t 20| . On the other hand, it has already been shown that the EPAC method works well in a one-dimensional 
^~«c double we,, s y s,e m |. This success suggests that the EPAC should he very useful for the calculation of 
quantum dynamics of the systems in which quantum coherence is significant. 

Most of actual potentials appearing in chemical systems and condensed matter systems have asymmetric sha pes . 
For example, asymmetric double well potentials have been widely used to describe the proton transfer reactions |2l| 
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or the first-order phase transitions [2^, [2^] . The other asymmetric potentials such as the Morse potential have also 
been very important and have been repeatedly used in chemical physics as a realistic model of the intramolecular 
vibration. 

In the present paper, the EPAC method is therefore applied to the evaluation of real time quantum correlation func- 
tions for one-dimensional asymmetric systems. We investigate whether the EPAC method can successfully reproduce 
the exact correlation functions for the Morse potential. We also show that the EPAC calculation for a certain type of 
asymmetric potentials can be reduced to that for the corresponding symmetric potentials. This property is called the 
decoupling property, which is quite useful for the discussion of the asymmetrization effect in the quan tum statistical 
systems, for instance, the enhancement or suppression of the transition rate in chemical reaction |2j|. Furthermore, 
the decoupling property can aid us in reducing the computational cost required on the PIMD/PIMC calculation. In 
this paper, we analytically prove that this property holds, followed by the numerical examination of it. 

In Sec. II, we summarize the definition of a couple of types of effective potentials and review our EPAC method. 
Moreover, we explicitly present the standard effective potential and the EPAC position autocorrelation function for 
a harmonic oscillator as the simplest example. The EPAC analyses of one-dimensional systems with asymmetric 
potentials are shown in Sec. III. After showing the decoupling property of linear terms appearing in the classical 
potentials, we present the EPAC results for the asymmetric harmonic/ anharmonic oscillators. The conclusions are 
given in Sec. IV. 



II. EFFECTIVE POTENTIAL ANALYTIC CONTINUATION METHOD 
A. Effective potentials and analytic continuation 

We hrst introduce the effective classical potential defined by Feynman [l|. Consider a quantum system where a 
quantum particle of mass m moves in a one-dimensional potential V(q) at temperature T. The quantum canonical 
partition function of this system is expressed in terms of the imaginary time path integral 

/•oo rqiPK)=q 

Zp = dq Vq e~ 8 *'\ (1) 
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where (3 — 1/fcgT and Se is the Euclidean action functional 
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S E [q] = dr 
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(2) 



-mq 2 + V(q) 



After a Fourier decomposition of the periodic paths q{r) = (1//3H) J^L-oo e luJ ™ T <i( UJ n) with the Matsubara frequen- 
cies ui n = 2Tm/ph, we get US 

1 poo 00 2 poo poo 

Z P = ^ dq Q n ^# d(Req MWmqM) e-°"\ (3) 



where qo is the path centroid (zero mode): qo = q(uJo)/(3h. Inserting an identity 1 = J^° aQ dq c S(qo — q c ) into the 
integral in Eq. J3Jl and integrating out all the Fourier modes, we obtain 
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dq c e-Vn^), (4) 

Here VZ(q c ) is the effective classical potential |l|, |9|, [10| , which is a function of the position centroid variable q c . Among 
the integration of all the Fourier modes, the zero mode integral dqo yields only the replacement qo — ► q c , so that 
the effective classical potential Vp(q c ) does not contain the fluctuation effect of the zero mode qo. Therefore the 
effective classical potential contains only the quantum fluctuation effects of q(uj n=j to) modes. 



Another type of effective potential [ill Il2| is briefly summarized next. At first, consider the quantum canonical 
partition function represented in terms of an imaginary time path integral in the presence of a constant external 
source J, 

/oo rq(0H)=q 
dq / Vq e~ SE l h e 0Jqo = e 0w ^ J \ (5) 

-oo •/ q(0)—q 

where w^J) is the generating function for the connected Green functions with zero energy. From Eqs. and J5|l, 
the generating function is written in terms of Vp (q c ), 



w p {J) = -log 



(0) 



/ m f°° da e-^^ e pJq - 
This can be rewritten in the phase-space centroid representation, 

(7) 
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dq c dpc e pj qa e -/3{^/2m+V|(g c )} 



where p c is the momentum centroid variable. The "conventional" effective potential Va(Q) is defined as the Legendre 
transform of the generating function wp{ J) 0,0], 

V P {Q) = sup {JQ~ Wfj(J) }. (8) 
j 

We refer to this effective potential as the standard effective potential. This is the leading order of the derivative 
expansion of the effective action, and contains the effects of both quantum fluctuation and thermal fluctuation in 
the quantum statistical system. The static properties of the system can be reproduced from Vp{Q) straightforwardly 
. For example, one can evaluate the expectation value of the position operator q using the relation 

(00 = Qmin, (9) 



trie qua 

HQ 



where <5 m i n is the position of the standard effective potential minimum. The value of Q m in is determined by the 
stationary condition dV/3/dQ\g = Q = 0. In the low temperature limit, the two effective potentials are identical 
while the ground state energy of the quantum system is given by the minimum value of Vp{Q) 0,0, 

E Q = (0\H\0) = lim V f3 (Q min ). (10) 

0^00 



On the other hand, one can also calculate the dynamical properties such as real ti me q uantum correlation functions 
from the standard effective potential Vp{Q). Using the effective action formalism [ll|,ll2|, at first imaginary time 
quantum correlation functions are directly obtained from the standard effective potential Vp{Q) with the leading order 
derivative expansion. Then by means of the analytic continuation [2^ real time quantum correlation functions can be 
readily obtained from the imaginary time quantities. This is the EPAC an approximation method based on the 
effective action formalism. Following the EPAC method, for example, the real time position autocorrelation function 
Cp(t) — (q(t)q(0))/3 can be expressed as Q 

C (t)^Cf(t) = (^^coth^jcos^-^^-^sin^ + g^, (11) 

where ujp is the effective frequency, 

V ^ Q=Q min 



4 



D) 
i_ 

CD 

c 




FIG. 1: The inverse temperature (5 dependence of Vp(q c — 0), Vp(Q = 0), and top for the quantum harmonic oscillator 
with the parameters h = fcs = m = u = 1. Note that of the harmonic oscillator is independent of the temperature. 

B. Example: harmonic oscillator 



To illustrate the characteristics of the standard effective potential Vr(Q) and the EPAC correlation function, 
here we explicitly present the results for a simple system: a quantum harmonic oscillator whose classical potential is 
given by 



In this system the effective classical potential can be written as 



1 22 1, fsmh(/3hw/2) 



Then we can evaluate the generating function wp(J) from Eqs. @ and (|14|l 

1 



w (J) 



2mu> 2 



t2 1, / . Bhw 
- -log (^2sinh^- 



to obtain eventually the standard effective potential as the Legendre transform of wpiyJ) 

Vp(Q) = ^muj 2 Q 2 + i log ^2 sinh . 



(13) 



(14) 



(15) 



(16) 



The minimum of Vp(ff) is located at the point Q = Q m in = 0, and the effective frequency is temperature- independent: 
LUp = to. This potential minimum Vp(Q — 0) is equal to the free energy of the quantum harmonic oscillator. 

In fact, in the low temperature limit, the potential minimum Vp(Q — 0) is equal to that of the effective classical 
potential VS(q c = 0), and they give the ground state energy of the quantum harmonic oscillator, 



hm V P {Q = 0) = lim V$(q c = 0) = 



(17) 



On the other hand, in the high temperature limit (/3 — > 0), these effective potentials exhibit quite different features. 
The standard effective potential minimum Vp(Q = 0) diverges to negative infinity, while the effective classical potential 
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minimum VS(q c = 0) converges to zero. This is because the standard effective potential Vp(Q) includes the effects of 
both quantum fluctuation and thermal fluctuation, while the effective classical potential VS(q c ) includes the effect of 
quantum fluctuation only. In Fig. \I]we show the inverse temperature /3 dependence of Vp(q c = 0), Vp(Q — 0), and 
up with the parameters h = kB=Tn = u> = l. 

Finally we obtain the EPAC position autocorrelation function for the quantum harmonic oscillator 

C% c (t) = ( 7-^— coth^^ ) cosujt - i ( — — | sin cut. (18) 

This i s eq ual to the exact correlation function. The EPAC method is thus exact for harmonic systems as well as the 
CMD 0. 



III. APPLICATION TO ASYMMETRIC SYSTEMS 



In our previous paper [8j, the usefulness of the EPAC method has been shown for a ^(g o — q) symmetric 
double well potential system in comparison with the CMD result. We expect that the EPAC method may also be 
useful for general anharmonic quantum systems. In the present section, we embark on a project to investigate how 
the EPAC method works in a Zi asymmetric system with the classical potential 0, , 

where natural units h = ks = m = 1 are employed. This potential is given by a polynomial approximation to a Morse 
potential, 

V(q) = § (l-e^f (20) 

regarding q as the small deviation from the minimum q = [3(3] . This polynomial approximation is very good for 
(3 > 10. The Morse pote ntial (|20|) corresponds to, for example, the parameters denoting the intramolecular vibration 
of the HC1 molecule [3JJ. For this system the inverse temperature /3 = 10 corresponds to the temperature T ~ 400 K. 
On the other hand, by a variable shift x = q + | , we can transform the potential (|19fl to the form 

V(x) = — --x + -x 2 + —x\ (21) 
v ; 64 4 8 100 ' v ' 

which is a potential asymmetrized by the linear term — | x only. Such additional linear terms exhibit interesting feature 
in the effective potential analyses; they are irrelevant both to the path integral calculation © and to the Legendre 
transformation (JSJ) . This property ensures that the standard effective potential for the asymmetric potential (|21|l can 
readily be given by that for the corresponding symmetric potential without the linear term. This is the decoupling 
property, which has been discussed for other types of effective potentials in the context of the nonperturbative 
renormalization group [3, . In the following we give a proof of this property for the standard effective potential, 
and then apply it to the evaluation of the standard effective potentials for an asymmetric harmonic oscillator and the 
asymmetric anharmonic oscillator (|21|l . 



A. Decoupling property of the additional linear term 



Consider a quantum system described by a classical potential of the form 



V(q) = U(q) + fq, 



(22) 
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where U(q) is any classical potential and / is a constant. The effective classical potential can be written as 



V$(q c ) = U C p{q c ) + fq c , 



(23) 



where U^{q c ) is the effective classical potential for U(q). This is because in the Euclidean action Se the linear term 
consists of only the zero mode go, 

r 0H 



dr fq(r) = f(3hq , 



(24) 



and never generate or receive any quantum correction when we integrate out all the Fourier modes to evaluate the 
effective classical potential. 

Next, from Eqs. © and H23[l the generating function for V(q) is given as 



wp{J) 



2TT/3H 2 



dq c e^ u ^e^ J -^ 



(25) 



Supposing that is the generating function for U(q), an evident relation wp(J) = uri (J — /) holds. Furthermore, 
supposing J*(Q) is the solution of the equation Q = dwp/dJ and J{j(Q) is that of Q = duip /dJ, we get a relation 
J*(Q) — JjjiQ) + /• The standard effective potential for V(q) can therefore be written as 



Vp{Q) = {JQ-wp{J)}\j=,,, 



(Q) 



Ju(Q)Q + fQ- «iog 
p 

Up(Q) + fQ, 
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dgcC -/JD|( a .) e /3^W)9« 



(26) 



where Up(Q) is the standard effective potential for U(q). Equations (|23|l and (|26|l ensure that the functional form of 
the linear term remains unchanged through the path integration and the Legendre transformation, and that it never 
generates any quantum or thermal correction to U(q). That is, the linear term in the classical potential decouples 
from the effective potential calculations. 

Using Eq. 1)26(1 . we can directly obtain the standard effective potential Vp(Q) from Up(Q). This means that we 
have only to carry out the path integration and the Legendre transformation for U (q) . This decoupling property is 
especially useful when U(q) is a Z2{q «-> — q) symmetric potential, because we can simply discuss the asymmetrization 
effect as the additional linear term effect on Up(Q). Furthermore, the Zi symmetry reduces the computational cost 
required on the PIMD/PIMC calculation for the Zi symmetric U(q). This is because the required cite number in q c 
space is reduced by half by the Z 2 symmetry ■ 



B. Asymmetric harmonic oscillator 

As a simple example of the system with the additional linear term (1221 , we consider an asymmetric harmonic 
oscillator 

V(q) = ~ mu\ 2 + fq. (27) 
From Eqs. I|16() and (|26|l wc readily obtain the standard effective potential 

Vp(Q) = \muj 2 Q 2 + fQ + ^ log ^2 sinh . (28) 

The minimum of Vp(Q) is at the point Q = Q mm = —//mo; 2 , and the effective frequency is still temperature- 
independent: ujp = lu. The value Q m i n coincides with the exact position expectation value (q)p, while the minimum 
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value, 



2mw 2 

gives the exact ground state energy in the low temperature limit 



Vf3(Q min ) = + - log 2 sinh ^— , (29) 



f 2 So; 

So = Hm Vp(Q min ) = -^-a + "V ( 30 ) 



The EPAC position autocorrelation function is then obtained as 

= f^coth^cos^-if-A-'jsinwi + ^-j. (31) 
This present result is also equal to the exact correlation function. 

C. Asymmetric anharmonic oscillator 

In this section we attempt to apply the EPAC method to the asymmetric anharmonic system with the potential 
(|19|l . To obtain the parameters Q min and up needed for the EPAC correlation function now we propose the 

following computational schemes. 



Scheme (A). By carrying out the path integration and the Legendre transformation for the classical asymmetric 
potential (|19|l . we obtain the standard effective potential Vp(Q). We then fix the parameters Q min and up by 
minimizing Vp(Q). 

Scheme (B). We first carry out the path integration and the Legendre transformation for the symmetric part in the 
converted potential l|21(l . 

PW = f + ^ + I5o* 4 ' < 32 > 

to obtain the symmetric standard effective potential Up(X). Using Eq. Ij26(l with / = — 4, we can write the asymmetric 
standard effective potential as 

%{X) = Up(X)-~X. (33) 

We then obtain the minimum point A m i n and the effective frequency up by minimizing Vp(X). The parameters Q m i n 
and cup are obtained from the relations Q m in = X m [ n — | and top = up. 



In the following we perform the numerical calculation with scheme (A) to test the EPAC method itself. The 
numerical calculation with scheme (B) is also performed to check the decoupling property l|26(l . 

To evaluate the standard effective potentials Vp(Q) in scheme (A) and Up{X) in scheme (B) using Eqs. © and JBJ), 
we need to precompute the effective classical potentials Vp(q c ) and Up(x c ), respectively 33]. Employing the PIMD 
technique and the fitting procedure described in Ref. 0, we calculated Vp(q c ) and Up(x c ) at the inverse temperature 
(3 = 0.1,1,10, and 100. We then calculated the standard effective potential Vp(Q) and Up(X) by carrying out the 
numerical integration © and the numerical Legendre transformation JSJ) . As for scheme (A) , Fig. [2] shows the 
evaluated effective potentials for the asymmetric potential V(q) [Eq. (|19fl ]. As seen in Fig. Ufa) and[5{b), at the 
low temperature f3 = 100, the standard effective potential Vp(Q) has almost the same shape as the effective classical 
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FIG. 2: The inverse temperature j3 dependence of the effective potentials for the asymmetric classical potential 1191 : (a) 
the effective classical potential Vp(q c ) and (b) the standard effective potential Vp{Q). In this plot we set Vg(0) = and 
V p (Q nin ) =0. 

potential Vg(q c ). This is a general property of the effective potentials 0, [3, [2]]] . With regard to Scheme (B), the 
evaluated effective potentials for the symmetric part U(x) [Eq. l|H2|l] are shown in Fig. 01 It is also seen in this figure 
that Up(x c ) and Up(X) have more similar shape at lower temperature, as in Fig. |21 

Next we evaluated the parameters Q m i n and up. With scheme (A) we directly obtained them from the computed 
Vp(Q), while in scheme (B) they were obtained from the evaluated effective potential Up(X) using the relations 
Qmin = X m i n — I and uip = (Dp. Figure 01 shows Q mm and u>p calculated by means of scheme (A) and (B), and also 
shows the exact result 

where the energy eigenstates |n) and the energy eigenvalues E n were obtained by a numerical integration of the 
Schrodinger equation. In this figure, we also plotted the effective frequency of the symmetric effective potential 
Up(X): LJg — y/ d 2 Up(X) /dX 2 \ x=Q . It is found that at each temperature the quantities computed in scheme (B) 
are in very good agreement with the corresponding ones in scheme (A). This indicates that the relation l|26(l actually 
holds in the numerical PIMD calculation and the numerical Legendre transformation. We can also see that each Q m i n 
agrees very well with the corresponding exact value (q)p; the effective potential approach to such a static property 
works successfully in this asymmetric system. 

Here we give a comment on the asymmetrization effect. As seen in Fig. each effective frequency u)p(= uip) 
is larger than at each temperature, that is, the asymmetrization effect enhances the effective frequency of the 
system. This property can be simply understood as follows: For the weakly anharmonic system, the symmetric 
standard effective potential Up(X) can be well approximated by a power series expansion to the fourth order, Up{X) = 
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FIG. 3: The inverse temperature (3 dependence of the effective potentials for the symmetric classical potential (13211 : (a) the 
effective classical potential Up(x c ) and (b) the standard effective potential Up(X). In this plot we set £7(0) = 0, Up(0) — 0, 
and Up(X mia ) = 0. 
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FIG. 4: The inverse temperature (3 dependence of the minimum point Q m i n and the effective frequency u)g. The exact position 
expectation value (q)p and the effective frequency of the symmetric effective potential top are also plotted. 
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C + ^ujp 2 X 2 + jApX 4 with a constant C and a positive coupling constant Ap. Then the asymmetric standard effective 
potential can be written as Vp(X) = C — \X + \{jj s ^X 2 + jApX 4 . Supposing A m ; n is the minimum point of Vp(X), we 
obtain the effective frequency of Vp(X) as Up = \J^p 2 + 3ApX^ in . Since the parameters ujp 2 and Ap are positive, we 
get an inequality u)p > w|, where the equality holds in the harmonic limit Ap — > 0. Figure0]also shows that u>p{= dip) 
becomes closer to ujp as the temperature increases. This implies that the standard effective potentials Up(X), Vp(X), 
and Vp(Q) become effectively harmonic ones at the higher temperature. 

Finally we constructed the EPAC position autocorrelation functions C^ c (t) using the computed quantities Q m i n and 
uip. Figure shows the real part of C^ c (t) together with the real part of the exact quantum position autocorrelation 
function, 

^( < ) = ^EE e ^ £ " e " l(Bm " B " )t/;i IH^}| 2 , (35) 



at the inverse temperature j3 — 0.1, 1, and 10. We can see that each EPAC correlation function reproduces well the 
exact correlation function at t = 0; the static property Cp(0) can be approximated very well by only two quantities 
Qmin and cop, 



Cp(0) ~ C£ c (0) = ^- coth + Q min . (36) 



This means that the leading order derivative expansion employed in the EPAC 8] is very good in the calculation of 
the static property Cp(0). The value (7g C (0) can be improved by employing the higher order derivative expansion; 
it should become exact in the infinite order expansion if the expansion is a converging one. On the other hand, as 
the time t increases, each EPAC correlation function deviates from the exact one. Such deviation becomes larger at 
the higher temperature. In particular, at (3 = 0.1 [Fig. |Sfa)], the exact correlation function Cp(t) damps rapidly 
with time, whereas the EPAC correlation function C^ c (t) shows stable oscillation. This damping behavior comes 
from the fact that the exact correlation function l|35l) consists of the many oscillation modes with the frequencies 
^m.n — (E m — E n )/h, and oscillation modes with high energies contribute to Eq. (|35[l more at higher temperature. 
On the other hand, the EPAC correlation function consists of only one oscillation mode with the frequency wp at any 
temperature, so that it always exhibits a single-mode oscillation. As the temperature lowers [Fig. Efb) and Etc)], the 
EPAC correlation function C^ c (t) becomes closer to the exact Cp(t). This is because at lower temperature, almost 
only one oscillation mode contributes to Eq. I|35(l to make the exact Cp(t) show a quasi-single- mode oscillation. This 
property generally holds for a weakly anharmonic system 

BBS 

though it does not hold for strongly anharmonic 

systems such as deep double well systems js^ . 

The results shown in Fig. [S] imply that the leading order derivative expansion employed in the EPAC method 
becomes worse at longer time and at higher temperature. To reproduce the damping behavior of the exact Cp(t) at 
high temperature, we need the higher order derivative expansion of the effective action. However, our calculation at 
lower temperature (3 — 10 [Fig. [5[c)] corresponds to, for example, the dynamics of the intramolecular vibration of the 
HC1 molecule at temperature T ~ 400 K. Therefore our EPAC method even in the present form (|llfl should work 
very well for such chemical systems in the realistic temperature region. 

Finally we mention the applicability of the EPAC method to a nonpolynomial potential, in particular, an asymptot- 
ically sublinear potential such as lim^^^ V(q)/\q\ < 0. Most of the actual potentials appearing in chemical systems 
and condensed matter systems have this type of potential shape, for example, the Morse potential, the Lenard- Jones 
potential, and other intramolecular and intermolecular potentials. An important property of these potentials is that 
they represent the dissociation in the long distance limit, ft is quite nontrivial to make the dissociation effect involved 
in the effective potential approach. In fact, some difficulties have already been pointed out 01 • From the chemical 
interest, this is a challenging subject to be examined. 
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FIG. 5: The real part of the exact quantum position autocorrelation function Cp{t) and the real part of the EPAC position 
autocorrelation function Cp C (t): (a) at (3 = 0.1, (b) at (3 = 1, and (c) at (3 = 10. 

IV. CONCLUDING REMARKS 

In the present work we have focused on the EPAC analyses of one-dimensional asymmetric systems. At first, for 
the asymmetric harmonic oscillator (|27ll , we have analytically shown that the exact real time position autocorrelation 
function is derived by means of the EPAC method. As for the asymmetric anharmonic oscillator Ijl9(l . we could obtain 
the EPAC correlation function numerically. In this case, at first we have computed the standard effective potential 
Vp(Q) by means of the PIMD technique and the numerical Legendre transformation. Then we have obtained the 
position of the effective potential minimum Qmin an d the effective frequency top which are needed to construct the 
EPAC position autocorrelation function C^ c (t). It has been shown that the static properties of this asymmetric 
system, (q)p and Cfj(0)(— (q(0)q(0))p), are well reproduced by Q m i n and C^ c (0), respectively. As for a long time 
behavior, we have found that at lower temperature the EPAC correlation function Cp C {t) agrees with the exact Cp{t) 
very well, whereas the EPAC result becomes worse as the temperature increases. This is because the EPAC method 
approximates the multimode contribution to the exact real time correlation function by the single-mode contribution 
with the effective frequency top. In a chemical aspect, however, it is quite remarkable that a good agreement of the 
EPAC correlation function with the exact one is achieved for the realistic potential in the realistic temperature region. 

We have also shown that a linear term in a classical potential is irrelevant to the effective potential calculation. For 
a classical potential V(q) — U{q) + fq, the standard effective potential is given by Vp(Q) = Up(Q) + fQ. We have 
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proved this decoupling property analytically, and have also checked it numerically by calculating Vp(Q) and Up(Q) 
independently. 

In the present work, it has been shown that the EPAC method is practically applicable to simple systems with single 
degree of freedom. More important applications in chemical physics are, however, to the systems with many degrees 
of freedom. For such applications, there are a couple of obstacles to be overcome. The first one is the computation of 
the generating function wp(J). To obtain wp(J) using Eq. @, we need to precompute the effective classical potential 
Vg(q c ) itself. However, this precomputation evidently demands more computational efforts as the system are more 
complicated. For many-body systems, Eq. JJJ) rather than Eq. © will be useful for more efficient evaluation of 
Wp(J), because the CMD sampling can be used to calculate the statistical-mechanical average of e l3Jqc appearing in 
the right-hand side of Eq. 33] . That is, we can use the numerical techniques developed in the CMD calculations 
for many-body systems UJi l2i| • The second obstacle is the evaluation of the standard effective potential Vp{Q) in 
many-body systems, i.e., the multidimensional Legendre transformation oiw@(J) [Eq. (JHJ)]. This transformation will 
also require heavier computations for larger systems. It is quite important to solve these numerical problems lying in 
the practical calculations for many-body systems. We are presently developing the numerical schemes along this line. 
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